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ABSTRACT 

A commonly accepted mechanism of generating baryon asymmetry in the minimal super- 
symmetric standard model (MSSM) depends on the CP violating relative phase between 
the gaugino mass and the Higgsino fi term. The direct constraint on this phase comes 
from the limit of electric dipole moments (EDM) of various light fermions. To avoid such 
a constraint, a scheme which assumes the first two generation sfermions are very heavy 
is usually evoked to suppress the one-loop EDM contributions. We point out that un- 
der such a scheme the most severe constraint may come from a new contribution to the 
electric dipole moments of the electron, the neutron or atoms via the chargino sector at 
the two-loop level. As a result, the allowed parameter space for baryogenesis in MSSM is 
severely constrained, independent of masses of the first two generation sfermions. 

PACS numbers: 11.30.Er,11.30.Fs, 12.60. Jv, 98.80. Cq 



While the Standard Model of particle physics continues to accurately describe a wide 
array of experimental tests many physicists suspect that the next generation of a unified 
field theory will be supersymmetric. This supersymmetric theory in its simplest form, 
MSSM[Q], may help to solve many of the outstanding problems in the Standard Model. 
Two examples of this sort are the coupling-constant-unification problem and the observed 
baryon asymmetry of the universe (BAU). It is the latter of these two that will be discussed 
in this paper. 

It has been demonstrated that SM is insufficient in generating large enough BAUj2|]. 
Particles lighter in mass but stronger in coupling are needed to make the electroweak 
transition more first order. Additionally, a new CP violating phase is required to generate 
enough BAU. It is very appealing that MSSM naturally provides a solution to both 
requirements 0. 

The top-quark partner, stop, which is naturally lighter than the other squarks can 
make the transition more first order, while there are plenty of new CP violating phases at 
our disposal in the soft SUSY breaking sector. In particular, it was shown that the most 
likely scenario is to make use of the relative phase between the soft SUSY breaking gaugino 
mass and the /i term of the Higgsino sector 0. In such a case, the BAU is generated 
through the scattering of the charginos with the bubble wall. The CP violation is provided 
by the chargino mixing. It turns out that in most parameter space of MSSM, a nearly 
maximal CP violating phase is needed to generate enough BAU. One immediate question 
is whether or not such a new source of CP violation is already severely experimentally 
constrained. It is not surprising that the most severe constraints are provided by the 
current experimental limits of the electric dipole moment (EDM) of the electron (d e ) and 
the neutron (d n ). 

Fortunately the lowest order (one-loop) contributions to various EDM's through the 
chargino mixing can be easily suppressed by demanding that the first two generations 
of sfermions to be heavier than the third oneQ, p|. For example, if one requires these 
sfermions to be heavier than 10 TeV, the one-loop induced EDM's will be safely small ||. 
In fact, such a scenario can even be generated naturally in a more basic scheme referred 
to as the more minimal SUSY model [0]. However, despite the enlarged parameter space 
of MSSM, thanks to all the intricate limits provided by accumulated data from various 
collider experiments, there is only a small region of parameters left within MSSM for such 
baryogenesis to workQ. 

In this letter we wish to point out that even if sfermions of the first two generations 
are assumed to be very heavy, there are important contributions to the EDM of the 
electron at the two-loop level via the chargino sector that strongly constraint the chargino 
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sector as the source for BAU in MSSM. Similar contributions to quark EDM also exist 
but the resulting constraint turns out to be relatively weaker. While this is not the 
first time that two-loop contributions are found to be more important than the one-loop 
ones,[§, |9], [K| [0], [12]], this chargino contribution and its relevance to BAU was never 
treated fully. 

In the case of chargino contributions, the two-loop contribution is dominant because 
the one-loop contribution is suppressed when the sfermions are heavy. This aspect is 
similar to those in Ref. ||, 0. In addition, the present case of the large CP violating phase 
in the chargino mixing and the light Higgs scalar, which is necessary to obtain a large 
baryon asymmetry, is also the same cause of a large EDM. Therefore, the resulting severe 
EDM constraint is very difficult to avoid in the mechanism of the chargino baryogensis 
by tuning parameters. 

The Model and Couplings 

Before we outline the physics of the chargino mixing in supersymmetric models we will 
set forth our conventions. We assume the minimal set of two Higgs doublets. Let the 
superfield $^ (Y = —1) couple to the <i-type field, $ u (Y = 1) to the w-type (see Ref. [|TT|1 
for our convention). The chargino fields are combinations of those of the wino (u^ #) and 
the higgsino (h+ L dR ). Denote ip L = (w£ , h+ L ) T , and %j) R = (w£ , K R ). The chargino 
mass terms, — = ipRM c ipL in our convention, becomes 



Where M 2 is the SUl(2) gaugino mass. Note that we choose a CP violating complex 
Higgsino mass fie 1 ^. The scalar components H u , H d of $ u , $^ have real vev's v u j v^2, v d / \[2 
respectively, and tan/3 = v u /vd- 

We use the bi-unitary transformation to obtain the diagonal mass matrix M D = 
U'M C U^ with eigenvalues m xl ,m X2 for the eigenfields xi-> X2- SUSY the The CP violating 
chargino mixing can contribute to the fermion EDM through the chargino-sfermion loop. 
Detailed analysis of such contributions can be found in the literature ||. As noted in the 
introduction, such contributions can be tuned to be small by making sfermions heavy|| 
(typically of 10 TeV or larger). Here we are interested in contributions to the EDM of a 
fermion that are still important even with very heavy sfermions. For this we find that the 
leading contribution is from diagrams of the type in Fig. 1. 

To evaluate the diagram, we exam gauge couplings of the Higgs bosons, H® = {y q + 
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\/2M w cos /3 




(1) 
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<P 9 )/V2, 

£y = JE*W + uMj^XjL + H.c. (2) 

ij 

Only the diagonal couplings in the chargino basis are relevant to the simple diagrams in 
Fig. 1 mediated by an internal photon. Therefore we define 

9i U = 9iu + Wfu = ^J^iJ^ih i 9i d = did + Wid = ^Pib^iu) ■ (3) 




/*(?) 



Fig. 1. A two- loop diagram of the EDM of the electron, or quarks. The chargino runs in 

the inner loop. 

The complex mixing amplitudes are written in terms of the real couplings g s and g p . 
In the same spirit, the complex neutral Higgs fields are decomposed into the real and 
imaginary components (p® — h G q + iaP (q = u,d). Note that h° d and mix in a CP 
conserving fashion at the tree- level, so are aP u and a d . 

*,)-'{%)■{'*)-{%)■ 

(cos a — sin a \ ( sin (3 — cos (3 \ 

' S = \ R ■ R " (5) 

sin a cos a J \ cos p sm p J 

The EDM calculation involves the Higgs boson propagators which are defined as 

(W^ = i E ^J/(P 2 " M l) , <^>p 2 = < E ^/(P 2 - M a) • (6) 
The Z factors can be shown to be real at the leading order with the explicit form, 

7 d,d _ 7 u,u _ 2 yd,d _ yu,u , o a 

L ±,H — L ±,h — cos a > Zj ±,G — Zj ±,A — ± cos P i 

Z ±,h = ±,h = sm a , Z ±,A = z ±,g = ± sm p , 

Z ±,H = 2 S111 2ft = _Z ±U ' Z ±U = ± 2 S111 2 ^ = _ ^±,G 

Z£"=2%i tor a = h,H,A,G . 
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For completeness, our list includes the unphysical Goldstone boson G°, which does not 
contribute to EDM. Other sum rules are 

E Z™J = 25™'5 S>+ . (7) 

a=hHAG 



The electron EDM via the Fig. 1 is given by 
I d e \ ot gvn e ^ g. 



e J 167T 3 M w cos (3 ^ q m Xi 



2 \ / 2 \ / 2 

mz . \ „ a I rnz . \ „ ^ / m: 



Here the Barr-Zee0 functions are defined as 

Kn(z) = Z -C y" lnV dy f(z) = Ro{z) _ 2Ki ( z) + 2K ^ z) g{z) = Kq ( z) (Q) 

2 Jo y(l -y)-z 

For the EDM of the down quark, we simply use the charge ratio | to give (dd/e) = 
^(d e /e)(m d /m e ). While for the EDM of the up quark, we need to replace Z q * d — > Z q,u 
in Eq. (||) as well as the obvious charge ratio — | and replacement of m e — > m u . In the 
Appendix, we offer a more compact analytic form of these results together with additional 
details which include the radiative correction to the Higgs mass in the simplified form 
suggested in Ref . JTjJ| . 



Since the charginos do not couple to the gluon, there is no chro mo-EDM generated |TT| 



Note that if one wishes to include the contribution with the internal photon replaced 
by the Z boson, it is necessary to include the off-diagonal chargino couplings of the Z 
and the Higgs bosons. We ignore such contributions here because they are expected to 
be much smaller than that of the photon which was confirmed in previous similar two 
loop calculations |jlO|. In particular, the electron EDM via Z is highly suppressed by the 



small value of the Z vectorial coupling to the electron due to the approximate relation 



sin 2 6y/ w j. There are other two loop diagrams with CP violation originated from the 
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same phase such as the ones with chargino-neutralino loop mediated r jH + W~ effective 
vertex or ~fW + W~ (W EDM) effective vertex. We do not include them here because these 
contributions are expected to be small (by roughly an order of magnitude) as suggested 



by previous two loop calculations [|ll], |T2[ . In any case, these additional diagrams form a 
separate gauge independent set. 

Because the imaginary parts of the off-diagonal entries in Mq are zero in our conven- 
tion, we obtained the following sum rules: 

E 9lu m x, = ~^{U' ] M D UU = , E 9i A m Xi = . (10) 

i i 

Therefore, g^ q = —gi, q ( m xi/ m x2)- ^ ^ s eas y ^° see *^ a ^ ^ n case °f degenerate masses 
m xi = m X2 , perfect cancellation occurs yielding zero EDM. 
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Based upon another fact that the diagonal scalar coupling of XiG°Xi is zero, we can 
show that sin (3gf u = cos f3gf d . Therefore, each of the four CP violating coefficients gf 
can be simply related to one of them, say g[ u , which again depends on the fundamental 
MSSM parameters, tan /3, fie^, M\, M 2 . The usual SUSY breaking terms include the last 
two parameters as well as the trilinear-sfermion-coupling, the A term, which is not relevant 
in the our analysis because it does not participate directly in this particular mechanism 
of baryogenesis||. If we replace charginos by stops in the inner-loop, the effect of the 
relative phase of A and /i can contribute to the two-loop EDM as studied in Ref.||. The 
stop-loop effect can be small if A t is small, if A t is in phase with /x, or if the left-handed 
stop is very heavy but the right-handed stop is rather light. This last scenario is preferred 
by BAU. Such a large mass gap will suppress stop-mixing and kill the EDM contribution 
via the stop-loop. In addition, it has been concluded by many groups|3[] that using CP 
violating mixing of the stop to generate BAU is much more difficult than using that of 
the chargino. 



Numerical analysis and baryogenesis 

To our current knowledge, the experimental constraint on the electron EDM has become 
very restrictive: 

|4| < 1.6 x l(r 27 e cm (90% C.L. Ref.[13]) . (11) 

Since the tree- level Higgs mass relation |l| predicts a light Higgs m h o < m z that has 
already been ruled out by experimental searches at LEP 2, our analysis has included the 
leading mass correction |T3| at the one- loop level. For completeness, the resulting Higgs 
mass dependence on tan (3 in this scheme is illustrated in Fig. 2. Fig. 3 shows the tan /3 
dependence of the predicted value of the electron EDM from different contributions due 
to the Higgs bosons, A , H° and h°. We show the case of maximal CP violation when 
(f) = it/2, as required by baryogenesis [pT^j , with masses at the electroweak scale, Ma = 150 
GeV, M.2 = \i = 200 GeV. Note that, in this case, the h contribution dominates until about 
tan/3 m 3. The H contribution becomes dominant for tan/? > 5.4. When tan/3 becomes 
large, the increase of the Yukawa coupling of the electron overwhelms the reduction of CP 
violation in the chargino sector. This gives the rise of the electron EDM as tan (3 increases. 
The same effect happens to the EDM of the ci-quark, but not the u-quark. Fig. 4 shows 
the electron EDM contour plot versus M 2 and fi for the case tan f3 — 3, M4 = 100 GeV, 
and = 7r/2. In the many calculations of BAU in MSSMj3| the largest uncertainty seems 
to come from the calculation of the source term for the diffusion equations which couples 
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to the left-handed quarks jig, [L9|. Using the latest summary of the situation in Ref . P0j] as 
a reference point, large BAU (2 < r/ 10 = {n B — n^)/n^ x 10 10 < 3) requires tan (3 < 3 with 
the wall velocity and the wall width close to their optimal values v w ~ 0.02, ~ 6/T, 
fi ~ M2 and CP phase sin0 close to one. Note that a smaller tan/5 gives a larger BAU, 
however, it tends to give a small lightest Higgs mass which violates the LEP II limit unless 
the left stop is much heavier than 1 TeV. Using the SUSY parameters in the above range, 
the numerical analysis in our figures indicates that the predicted value of the electron 
EDM is more than a factor of 5 to 10 bigger than the experimental limit on the electron 
EDM in most of the BAU preferred parameter range. In fact, if sin</> = 1 and tan/? = 3, 
then the parameter space allowed by the electron EDM limit is limited to a narrow strip 
with fi ~ M2 and /1 has to be as large as 600 GeV in order to satisfy this EDM constraint. 
The range of values for /i and M2 (both smaller than 250 GeV) presented in Ref. |20 
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all ruled out. Unless the numerical constraint on BAU in Ref.|20[ is relaxed by an order 
of magnitude, it seems to be very difficult for the chargino mechanism for BAU to be 
compatible with the electron EDM constraint. 

On the the other hand, for the neutron EDM, our analysis indicates the current ex- 
perimental limit in Eq.(12) gives only marginal constraint on MSSM parameters required 
for chargino BAU. 

With the quark EDM, one uses the quark model to predict the neutron EDM. A new 
limit |L5| , \d n \ < 6.3 x 10~ 26 e cm (95% C.L.), for the neutron EDM has been reported 
based on the combination of the recent data of low statistical accuracy and the earlier 



measurement JT6|] . This combination of the old and the new results has been criticized 
in Ref. [T7J] . As shown in the contour plot of Fig. 5, using the parameters suggested by 
the chargino baryogenesis mechanism, our predicted EDM value is around the size of the 
more conservative experimental limit: \d n \ ~ 12 x 10~ 26 e cm, recommended in Ref.|17 . 



Due to large theoretical uncertainties in the relation between the quark EDM and the 
neutron EDM, the constraint from neutron EDM on the parameter space cannot be as 
important as that from the electron EDM even if the more stringent limit is used. 

Note, however, that the uncertainties in the calculation of non-equilibrium electroweak 
baryogenesis process is far from settle. For example, in the latest review by the group in 
Ref|H] a small CP violating phase of 10~ 2 may be sufficient to generate BAU. In that 
case even the larger value of tan f3 is allowed. For this purpose, in Fig. 6, we also plot the 
electron EDM for tan (3 up to 50. 
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Conclusion 



The baryogenesis in MSSM requires the lightest Higgs boson to be light in order to get 
a strong first order phase transition. It also requires the CP violating phase in chargino 
mixing to be large in order to get large enough BAU. As we have discussed, both re- 
quirements imply the predicted values of the EDM's of the electron and the neutron to 
be large. For sin0 = 1 and tan (3 = 3, the current electron EDM constraint requires 
fi ~ M 2 — 600 GeV. Taking the uncertainty in the calculations of BAU in the literature 
into account, it is probably still premature to claim that this particular mechanism of 
baryogenesis is absolutely ruled out, but it is clear that the precision measurements of the 
EDM of fermions, especially the electron EDM, give a tight constraint on the mechanism. 
Note added: While the paper is in referee process we receive a manuscript ( |hep 



ph/0207277| ) with calculations that overlap with ours. Our numerical results agree with 
this later calculation. 
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Appendix: Higgs Potential with radiative corrections 
in the MSSM and Electric Dipole Moments 

The Higgs potential has the form 

V = m 2 H jH d \ 2 + m 2 H jH u \ 2 + (-m 2 2 H d H u + H.c. ) 
+ M + 9l){\H d ?-\H u \ 2 ) 2 + r\H u \* + ... 

At the tree level, SUSY requires the dim=4 coefficient r = 0. However, it arises from the 
large top-stop-loop correction. Denote 

(H d ) = V d , (H U ) = V U , V 2 = V 2 + V 2 
tan/3 = V u /V d , m 2 w = \g 2 V 2 , m% = \{g 2 + g 2 2 )V 2 . 

We try to derive the mass matrix of the CP-even Higgs bosons, which correspond to 
the real part of the complex fields. We use superscripts R, I to abbreviate the real and 
imaginary parts. The first derivatives of the potential are 

{dV/dHf) = 2m 2 H H^-2m 2 2 H^ + l(g 2 + g 2 )(\H d \ 2 -\H u \ 2 )H^ 
(dV/dH*) = 2m% u H* - 2m&H? - \{g\ + gi)(\H d \ 2 - \H U \*)H* + 4r|F u | 3 ' 
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The minimization condition can then be written as 



tan (3 + |m| cos 2/5 = 



m# - m 2 2 cot (3 - ~m| cos 2/5 + 2r\^ 2 sin /3 = 



(15) 



Continue to obtain the second derivatives, 



{d 2 V/dHP) = 2m\ 2 tan (3 + 2M|c| 
d 2 V/{dH*dHf) = -2mj 2 -m| sin 2/5 , (16) 

(d 2 V/dH* 2 ) = 2m\ 2 cot /5 + 2s 2 (Mf + 4rV 2 ) 

(d 2 V/dH 12 ) = 2m 2 2 tan/5 
<9 2 V/(<9^<9#j) = 2m\ 2 . (17) 

{d 2 V/dH™) = 2m 2 2 cot/5 

The basis defined in Eqs.(f|||) agrees with that in Martin's review JTJ. One can easily show 
that G is massless as it is the unphysical Goldstone boson. The mass of the pseudoscalar 
A is 

m A o = 2m 12 / sin 2/5 , m H ± = m A o + m w . (18) 

The coefficient m\ 2 corresponds to the non-hermitean quadratic term in the Higgs poten- 
tial. If Tn 2 2 — 0, the Lagrangian possess a Peccei-Quinn symmetry and it quarantees that 
M A o =0. It is practical to express all other masses in terms of m A o. From the second 
derivatives above, the tree-level mass matrix of the scalar Higgs bosons in the basis of 
h^, h Q d becomes 

^ 2 / m 2 A0 cos 2 (3 + m| sin 2 (3 — (m A0 + m%) sin/5cos/5\ 
\ — (m 2 A o + m 2 z ) sin /3 cos /3 m 2 A0 sin 2 (3 + m 2 z cos 2 (3 J 

where the subscript indicates tree- level quantities. One can then prove that (rriho)o < 
mz\ cos 2/3 1 . 

The leading correction from top-stop loops is 

M 2 1LT ^M 2 + T 2 ( 1 ° ) , T 2 = 4rV 2 4 = 2 39 2 m < 2 In {m~ tL m~Jm 2 ) . (20) 
\ / o7H?7% sin /5 v ' 



This formula can be found in Ref.|[14j|, where different schemes of approximation were 
studied. As we have uncertaintly from the SUSY breaking scale, it may be overboard to 
use the full-fledge 1-loop calculation. We use this leading approximation in the remaining 
study. The CP-even Higgs mass-squared eigenvalues are then given by 

1 



m ^o 5 /iO - 2 



M 2 n + M\ 2 ± J[M 2 U - M 2 22 } 2 + 4(M 2 ^ 2 



l l2) 



(21) 
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The mass of h° has been substantially raised above the tree level prediction which is lower 
than the experimental constraint. The corresponding mixing angle a is given by 



sin 2a 



cos 2a 



2M 2 



12 



[M 2 



ii 



MP 2 



l 22j 



4(M 2 -^ 2 



L 12) 



M 2 22 - M 2 



ii 



[M 2 n -M 2 22 ] 2 + 4(M 2 12 ) 2 



The eigen-masses (m 2 H0 > m 2 h0 ) are given by 



[22] 



m 2 H0 + m 2 h0 



\m 



//o - rn 2 h o) 2 



m 2 A0 +m 2 z + T 2 , 

[{m 2 A0 - m 2 z ) cos 2f3 + T 2 } 2 + (mi + m|) 2 sin 2 2/3 



In terms of these masses, the mixing angle a is determined at tree-level by 



sin 2a 



m 2 A0 + m| 

2 2~ ' 

Wl^O m h° 



cos 2a 



(m| -m 2 Jcos2/3-T 2 
2" 2" 

771 j^o ^feO 



(23) 



(24) 



sin 2(3 

From the vanishing of diagonal scalar coupling of we have spgf u = cpg( d for each 

mass eigenstate i. Therefore 



E<? g^qZ+fi = 9fA Z + d h + tan P Z ^ d h ) = gf )U {-\ sin 2a + tan (5 sin 2 a) 
= k>tan^[l - (m 2 - 4c 2 m 2 - m| - T 2 )/ (m 2 H - m 2 )] 



(25) 



and 



E g gf, q Z q + tH = g^ u ( z l\ H + tan = ^ u (| sin 2a + tan /3 cos 2 a) 
= k>tan/?[l + (m 2 A - 4c 2 m 2 A - m| - T 2 )/^ 2 , - m 2 )] 



I] g^qZ+^A = g?, u ( z + d A + tan Z^+i) = sin 2 /5 + tan I 3 sin2 = tan ^ 



(26) 



(27) 



The 2-loop EDM of the electron with the leading 1-loop mass correction becomes 



tL 



+ 1 



a 



gm e 



16tt 3 2M w cos /3 



gf u m Xl tan/3 



i | T 2 + M 2 +M 2 A {l + 2c w ) \ g(m 2 jM 2 ) 



m 



H 



mi 



rn 



T 2 + M 2 +M 2 (l + 2c 2/3 ) \ g(m 2 Xl /M 2 H ) + J{m 2 jM\) 



m 



H 



mt 



in 



XI 



in 



XI 



- (m Xl -> m X2 ) 



(2? 
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Fig. 2 The mass of the light Higgs boson h° versus tan (3. The lower set of curves corre- 
sponds to the tree- level result. The upper set of curves includes the leading one-loop 
(t,t) effect, for m^ L = 1 TeV and m~ tR = l^QGeV . Curves within each set are in the 
order of cases wia = 150, 200, 250, 300 GeV, from bottom to top. 
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Fig. 3 The predicted value of the electron EDM versus tan (3 from different contribtuions 
due to the Higgs bosons h°, A and H°, at the maximal CP violation when = it/2. 
Masses are set at the electroweak scale, Ma = 150 GeV, M 2 — ji — 200 GeV. 
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Fig. 6 The predicted value of the electron EDM versus large tan (3, at the maximal CP 
violation when = n/2. Masses are set at the electroweak scale, M 2 — fx — 200 
GeV. Curves from top to bottom are in the order of cases tua = 150, 300, 450, 600 
GeV. 
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